Hydro dynamic elliptic and triangular flow in Pb— Pb collisions at \/s = 2.76ATeV 
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It is shown that a simuhaneous comparison of both eUiptic and triangular flow from (2+1)- 
dimensional viscous fluid dynamics with recent measurements in Pb+Pb colhsions at the Large 
Hadron CoUider (LHC) favors a smaU specific shear viscosity (t?/s)qgp ~ l/(47r) for the quark- 
gluon plasma. Using this viscosity value, the relative magnitude of the elliptic and triangular flow 
is well described with Monte-Carlo Glauber (MC-Glauber) initial conditions while Monte-Carlo 
Kharzeev-Levin-Nardi (MC-KLN) initial conditions require twice as large viscosity to reproduce the 
elliptic flow and then underpredict triangular flow by about 30%. 

PACS numbers: 25.75.-q, 12.38.Mh, 25.75.Ld, 24.10.Nz 



1. Introduction. Much attention has been given re- 
cently to the extraction of the shear viscosity to entropy 
density ratio (i.e. the specific shear viscosity rj/s) of the 
quark-gluon plasma (QGP) from elliptic flow data in rela- 
tivistic heavy-ion collisions A major road block in 

this effort is insufficient knowledge of the initial shape of 
the thcrmalizcd fireball created in these collisions, whose 
initial ellipticity is uncertain by about 20% jl7H2lj . This 
induces an 0(100%) uncertainty in the value of (??/s)qgp 
extracted from elliptic flow [J 0] . After the discovery 
of triangular flow in heavy ion collisions at Relativistic 
Heavy Ion Collider (RH IC) l22l - |23 | and Large Hadron 
Collider (LHC) energies j25l427| . followed by the confir- 
mation of its collective hydrodynamic nature [22, 28- 33j 
and the realization that shear viscosity suppresses higher 
order harmonic flow coefficients more strongly than el- 
liptic flow [l2l . \13L \29L |34| . ISSf , it was recently suggested 
[23l |25| . I36l438| that a combined analysis of the elliptic and 
triangular flow coefficients V2 and 173 could yield a more 
precise value for the QGP shear viscosity and thereby 
reduce or eliminate the model uncertainty in the initial 
deformation of the QGP fireball and its event- by- event 
fluctuations. This Letter presents such an analysis. 

Our study is based on a (2+l)-dimcnsional viscous 
hydrodynamic model with longitudinal boost-invariance, 
describing numerically the transverse evolution of the 
heavy-ion collision fireball near midrapidity. As in past 
work 0-S 3-211 we explore two different types of fluc- 
tuating initial conditions for the entropy and energy den- 
sity profiles, generated from Monte Carlo versions of the 
Glauber and KLN models (see [l^ and references therein 
for details of the implementation used here). 

The MC-KLN calculations were done using a Monte- 
Carlo sample of initial state profiles with identical prop- 
erties as those used in [3| for the calculation of trans- 
verse momentum spectra and elliptic flow in 2.76 A TeV 
Pb-Pb collisions at the LHC. For the x dependence of 
the gluon structure function in the MC-KLN model we 
used the power A = 0.28 [l^; the normalization factor 
for the initial entropy density was fixed by hand to re- 
produce the measured charged hadron multiplicity den- 



sity dNch/dr] for the 5% most central collisions [331; the 
measured dependence of dNct/drj on collision centrality 
(40| is then automatically reproduced reasonably well by 
the model [l^ (see Fig.[TJa)). MC-KLN runs were done 
with rj/s = 0.2 which, for this type of initial conditions, 
was shown to yield a good overall description of the mea- 
sured transverse momentum spectra and elliptic flow in 
200AGcV Au-Au collisions at RHIC ^ and gave an 
impressively accurate prediction for the unidentified and 
identified charged hadron spectra and elliptic flows in 



M 41^. 



2.76 A TcV Pb-Pb collisions at the LHC 

For the MC-Glauber runs we generated a new set 
of initial configurations that differ from those used for 
200AGcV Au-Au collisions in Q by the wounded nu- 
cleon to binary collision ratio. Taking the initial en- 
tropy density s{r±;b) = K( ^ n^^ {r±;b)+ xn^^ {r±;b)), 
we determine k and 2; by a two-parameter fit to the 
ALICE data shown in Fig. [Ija). Due to viscous 
entropy production during the hydrodynamic evolution, 
which itself depends on collision centrality, the fit value 
for X depends on the assumed shear viscosity. For MC- 
Glauber initial conditions we took rj/s = 0.08 since this 



value was shown in |13l . l25l . |35| to provide a reason 



able description of the charged hadron V2{pt) and v^^pr) 
data measured by the ALICE experiment; this results 
in a- = 0.118 for Pb-Pb collisions at the LHC. Both the 
MC-Glauber and MC-KLN initial conditions arc hydro- 
dynamically evolved with equation of state (EOS) s95p- 
PCE [m which matches numerical results from lattice 
QCD at high temperatures to a hadron resonance gas at 
low temperatures j43j and implements chemical freeze- 
out at Tchom = 165 MeV. The hydrodynamic output is 
converted to final hadron distributions along an isother- 
mal decoupling surface of temperature Tdcc = 120 MeV, 
using the Cooper-Frye prescription. 



In [21[ we showed that, due to similar fluctuation mech- 
anisms, the MC-KLN and MC-Glauber models generate 
similar third-order eccentricities £3 whereas the elliptic- 
ity £2, which is mostly controlled by collision geometry, 
is about 20% larger in the MC-KLN model. Event-by- 
event ideal 21 1 and viscous 38, hydrodynamic sim- 
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FIG. 1: (a) (Color online) Centrality dependence of charged particle multiplicity density as a function of A'part from the MC- 
Glauber (dashed) and MC-KLN (solid) models, compared with ALICE measurements [4Q] for 2.76 A TeV Pb-Pb collisions, (b) 
Charged particle pT-spectra from the MC-Glauber and MC-KLN models for different centralities. The most central (0—5%) 
results are compared with ALICE data [47l |. 



ulations with both realistically fluctuating l2 l |. | 4 2ll and 
doubly-deformed Gaussian initial conditions [38[ (with 
simultaneously non-zero £2 and £3 eccentricities) have 
shown that the hydrodynamic conversion efficiencies for 
translating initia l sp atial eccentricities into final flow 
anisotropies [s^, 13 ISIj although different for V2/e2 
and v^/e^, are very similar in the MC-KLN and MC- 
Glauber models. The similarities in £3 and differences in 
£2 between these models should thus straightforwardly 
reflect themselves in analogous differences in and ^3 
37 , 33 , allowing for an experimental distinction between 



the models. 

Evcnt-by-cvcnt viscous hydrodynamic simulations 
with full inclusion of unstable resonance decays arc at 
present numerically too costly for systematic flow studies 
over a range of viscosities, collision energies, centralities, 
and coUision systems. A recent event-by-event study by 
Schenke et al. [3^ with a restricted set of resonances 
showed that, compared to a full calculation, V2 (^3) was 
overprcdicted by 10-15% (25-30%). This is larger than 
the difference in these observables seen [21I 38, |42| be- 
tween cvent-by-cvent calculations and "single-shot hy- 
drodynamics" where the fluctuating initial conditions are 
averaged (after rotating each event by its second or third 
order participant plane angle [21I 38 1 to align the re- 
spective eccentricities) before the hydrodynamic evolu- 
tion instead of afterwards. For this reason we here use 
the single-shot approach, but include the full cascade of 
resonance decays in the flnal state. Our approach here 
differs from that in [2^ by replacing the singly-deformed 
Gaussian parametrization of the initial density used there 
by the ensemble-average of realistically fluctuating, non 



Gaussian initial proflles, and from [28j and [ij, [13|, |35| 



important effects of chemical non-equilibrium hadronic 
evolution on the elliptic flow 112 [46|. In 25| it was 
shown that, with the approach used in MC-KLN 
initial conditions with 77/s = 0.16 cannot describe the pt- 
integrated W3 measured in 2.76 A TeV Pb-Pb colhsions, 
whereas the MC-Glauber based event- by- event calcula- 
tions (with 77/3 = 0.08) of Schenke et al. [l3| appear to 
describe V3,(pt) at selected centralities reasonably well. 

2. Transverse momentum spectra. Figure[TJb) shows 
the charged hadron py-spectra for 2.76 A TeV Pb-Pb 
collisions at different centralities, for both MC-Glauber 
(77/5 = 0.08) and MC-KLN (77/5 = 0.2) initial conditions. 
For the most central (0—5%) collisions the spectra from 
both models agree well with published ALICE data. 
In more peripheral collisions the MC-KLN spectra are 
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FIG. 2: V2 and W3 vs. centrality, compared with ALICE U2{2}, 
«2{4}, W3{2}, and U3{4} data for 2.76^TeV Pb-fPb 
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FIG. 3: Eccentricity-scaled, pT-integrated 112,3 for the hydrodynamically evolved MC-KLN (a,b) and MC-GIauber (c,d) models, 
compared with ALICE V2,3 data for 2.76 A TeV Pb-Pb collisions [2^ scaled by their corresponding eccentricities (see text). 



harder than those from MC-Glauber initial conditions. 
This is a consequence of larger radial flow caused by 
larger transverse viscous pressure gradients in the MC- 
KLN case where the fluid is taken to have 2.5 times larger 
shear viscosity than for the MC-Glauber simulations, in 
order to obtain the same elliptic flow 0, Q . In periph- 
eral collisions these viscous effects are stronger than in 
more central collisions where the fireball is larger [3 ■ 
As shown in 
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49| , event- by- event evolution of fluctu- 
ating initial conditions generates, for small values of rj/s, 
flatter hadron spectra than single-shot hydrodynamics, 
especially in peripheral collisions, due to stronger radial 
flow driven by hot spots in the fluctuating initial states. 
Proper event-by-event evolution of the latter is there- 
fore expected to reduce the difference between the MC- 
Glauber and MC-KLN curves in Fig.[ljb) since this effect 
is relatively strong for 77/5 = 0.08 (MC-Glauber) [2l| but 
almost absent for r]/s = 0.2 (MC-KLN) 

3. pT-integrated elliptic and triangular Bow. In Fig- 
ure [2] we compare our pT-integrated V2 and V3 as func- 
tions of centrality with ALICE W2{2}, V2{4}, W3{2}, and 
■^214} data, extracted from 2- and 4-particle correlations 
[25|. For both models, V2,3 from averaged smooth ini- 
tial conditions lie between the experimental V2^3{2} and 
^2,3 {4} values. This is consistent with the theoretical ex- 
pectation 5^, m that Vn{2} (wti{4}) is shifted up (down) 
relative to the average flow by event-by-event flow fluc- 
tuations and was also found elsewhere 0, 0, • Upon 
closer inspection, however, and recalling that ideal single- 
shot hydrodynamics with smooth initial condition was 
shown [m to generate V2 similar to U2{2} from the cor- 
responding event- by- event evolution, it seems that the 
MC-KLN is favored since it produces V2 results closer 
to the 'y2{2} data. Unfortunately, a similar argument 
using V3 can be held against the MC-KLN model. To 
eliminate the interpretation difficulties associated with a 
comparison of average flows from single-shot evolution of 
averaged initial conditions with data affected irreducibly 



by naturally existing event-by-event fluctuations, we pro- 
ceed to a comparison of eccentricity-scaled flow coeffi- 
cients. 

Assuming linear response of U2.3 to their respective ec- 
centricities £2,3 (which was found to hold with reason- 
able accuracy for and W3 but not for higher order 
anisotropic flows [llj), we follow [lH and scale the flow 
W2,3 from single-shot hydrodynamics by the eccentricity 
£2.3 of the ensemble-averaged smooth initial energy den- 
sity, while scaling the experimental W2,3{2} and U2.3{4} 
data by the corresponding fluctuating eccentricity mea- 
sures e2,3{2} and e2,3{4}, respectively, calculated from 
the corresponding models. In [42| we justify this proce- 
dure for 1^2. 3{2} and W2{4} and also show that it fails for 
-(;3{4}/£3{4} since this ratio is found to differ strongly 
from v^/e^. 

The eccentricity-scaled elliptic and triangular flow co- 
efficients for the MC-KLN and MC-Glauber models are 
shown in Figs. [3lja,b) and[3ljc,d), respectively, and com- 
pared with the corresponding data from ALICE. The 
first thing to note is the impressively accurate agreement 
between the experimentally measured W2{2}/e2{2} and 
'y2{4}/e2{4}, showing that for elliptic fiow the idea of 
scaling "each flow with its own eccentricity" [s^l works 
very well. The same is not true for -(;3{2}/£3{2} and 
-(;3{4}/e3{4} for which the experimental do not at all 
agree (not shown), nor are they expected to ji^]. Sec- 
ondly, both W2{2}/e2{2} and W2{4}/£2{4} measured by 
ALICE agree well with the viscous hydrodynamic calcu- 
lations, for both the MC-Glauber and MC-KLN models, 
confirming that for each model the correct value oi rj/s 
has been used as far as elliptic flow is concerned. 

The bottom panels in Fig. [3] show the triangular flow 
V3. Clearly, with the viscosities needed to reproduce 
V2, the MC-KLN model badly disagrees with the ex- 
perimental data. The measured triangular flow is too 
big to accommodate a specific shear viscosity as large as 
0.2. Within the present approach, the only possibility to 
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FIG. 4: Eccentricity-scaled, pT-differential elliptic and triangular flow for 2.76 A TeV Pb-Pb collisions from viscous hydro- 
dynamics with MC-KLN (a,b) and MC-Glauber (c,d) initial conditions. The ALICE data [2^ are scaled according to their 
corresponding eccentricities, see text. 



avoid this conclusion is that somehow the MC-Glauber 
and MC-KLN models both underpredict the initial third- 
order eccentricity £3 by about 50%. With MC-Glauber 
initial conditions and 77/3 = 0.08, on the other hand, the 
ALICE data agree well with viscous hydrodynamics, even 
if the measured centrality dependence of v^{2} / e^{2} is 
slightly steeper than the calculated one. 

Summarizing Fig. [3l the ALICE data for the pr- 
integrated elliptic and triangular data strongly favor MC- 
Glauber initial conditions and, by implication, a small 
value of 77/5 ~ 0.08 for the specific QGP shear viscosity. 

4. pT-diffcrcntial elliptic and triangular How. Wc close 
this Letter by cross-checking, at one collision centrality 
(30—40%) where v^^pr) data are available |25[, the pr- 
differential anisotropic flows. The corresponding com- 
parison between data and theory is shown in Fig. |4l as in 
Fig. |3]we compare the eccentricity-scaled flows, plotting 
?^2,3/e2,3 for the models and V2{4}/£2{4} (v3{2}/e3{2}) 
for the elliptic (triangular) flow data. As seen in the 
upper panels, both initial state models describe the mea- 
sured elliptic flow well up to pr ^ 1 — 1.5 GcV/c; at larger 
Pt^ they overpredict vyApr) for charged particles - a 



problem noticed before [IJ, |16[ and possibly related to an 
imperfect model description of the measured final chem- 
ical composition [4l| . The disagreement at larger pT is 



worse for MC-Glauber initial conditions; this is likely re- 
lated to our earlier observation in Fig.[TJb) that our MC- 
Glauber pT-spectra are steeper than the MC-KLN ones 
in peripheral collisions - an artifact of our single-shot ap- 
proach and possibly remedied by a proper event- by- event 
hydrodynamical simulation. 

Figure Sl^b) shows again the disagreement between the- 
ory and experiment for triangular flow when we use MC- 
KLN initial conditions: the model strongly underpredicts 
the data at &\\pT, i.e. it gives the wrong slope for v^[pt)- 
With MC-Glauber initial conditions and correspondingly 
lower shear viscosity 77/s = 0.08 (Fig.UJ^d)), the measured 



V'iipr) is well described up to pt 1 GeV/c but overpre- 
dicted at larger pT- Again, the latter can be at least 
partially attributed to the fact that MC-Glaubcr pT- 
spectrum from our single-shot hydrodynamic approach is 
too steep at this collision centrality, which can in future 
studies be corrected by performing the hydrodynamic 
evolution properly event by event. 

5. Summary. Using a well-calibrated single-shot vis- 
cous hydrodynamic approach without hadronic cascade 
afterburner but properly implementing hadronic chem- 
ical freeze-out at Tchom ~ 165 McV and including a full 
set of resonance decays, we have shown that a combined 
analysis of the ALICE data for elliptic and triangular flow 
from 2.76 A TeV Pb-Pb collisions leads to a strong prefer- 
ence for initial conditions from the Monte-Carlo Glauber 
model, combined with a low value for the QGP shear vis- 
cosity 77/s~0.08, and disfavors the considerably larger 
viscosities of r]/ s ^ 0.2 that are required to reproduce the 
measured elliptic flow when assuming the more eccen- 
tric Monte-Carlo KLN initial profiles. Final confirma- 
tion of these conclusions will require a proper cvcnt-by- 
evcnt evolution of the fluctuating initial density profiles 
and coupling viscous hydrodynamics to a microscopic de- 
scription of the dilute late hadronic stage where viscous 
hydrodynamics breaks down [H^, and a similar analysis 
of recently published PHENIX data at lower RHIC ener- 
gies [H] ■ Given the large magnitude of the underpredic- 
tion W3 in the MC-KLN model observed here we doubt, 
however, that such more sophisticated approaches will be 
able to reverse the conclusions drawn here. 
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